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We compute the two-loop corrections of 0{a s ) to the Yukawa couplings in the 
framework of the Minimal Supersymmetric Standard Model (MSSM). The calcula- 
tion is performed using the effective Lagrangian approach under the approximation 
of neglecting the Higgs boson mass with respect to the top quark, gluino and all 
squark flavour masses. As an application we derive the 0(a 2 ) corrections to the 
partial decay width of the lightest Higgs boson to a bottom quark pair. We find 
that the two- loop corrections are sizable for large values of tan (3 and low CP-odd 
Higgs boson mass. With our calculation of the 0(a 2 ) corrections the remaining 
theoretical uncertainties reduce below a few percent. 
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1 Introduction 



One of the main purposes of the CERN Large Hadron Collider (LHC) is the search for 
Higgs bosons. The discovery of a light Higgs boson is a decisive test for all models 
predicting supersymmetric (SUSY) particles at the TeV scale and in particular of the 
MSSM pp. A remarkable feature of the MSSM is the restricted Higgs sector. This allows 
Higgs searches without any assumption about the mechanism of SUSY breaking, but only 
constraints from the Higgs sector [2j. For Higgs boson searches at the hadron colliders, 
two new complementary benchmark scenarios, the "small a e ff" and the "gluophobic" 
scenarios, have been proposed in addition to those used at the CERN Large Electron- 
Positron Collider (LEP) for the MSSM Higgs searches at hadron colliders. 
More precisely, in the "gluophobic" scenario the gluon fusion process is strongly suppressed 
due to cancellation between top quark and squark loop contributions. Nevertheless, the 
channel ti — > tih — > tibb is enhanced as compared to the Standard Model (SM) case, so 
that this becomes the most promising detection mode. In the "small a e g" scenario the 
decay width for h — > bb is much smaller than its SM value. In this case the complementary 
channel h — > 77 is enhanced as compared to the SM and it becomes the preferred detection 



mode. 

For a light Higgs boson (rrih < 130 GeV) the decay h — > bb is the dominant mode, 
but its detection at hadron colliders is difficult due to large QCD backgrounds for the b 
jets. However, at lepton colliders the Higgs boson search relies on b tagging that can be 
performed with high efficiency. 

For the discovery of the Higgs bosons, the cross section of the main production chan- 
nels, decay widths and branching ratios are necessary to be known with high accuracy. 
Within the SM the radiative corrections to the fermionic Higgs decay were intensively 
studied in the literature. The QCD and EW corrections are known up to the three-loop 
order: 0(a^) originating from the light degrees of freedom were first derived in Ref. [3] 
and the top-induced 0{a 3 s ) corrections in Ref. [I]; the QCD-EW interference contributions 
of 0(a^x t jj can be found in Ref. [5]. In this paper we concentrate on the fermionic Higgs 
decay in the MSSM, for moderate Higgs boson masses M h < 130 GeV. The processes 
h — >■ bb and h — > t + t— (the second most important decay mode) are affected by large 
radiative corrections for scenarios with large values of tan /3 and moderate values of the 
neutral CP-odd Higgs boson mass Ma (tan/3 > 20 , Ma < 250 GeV) [6]. Apart from pure 
QCD and EW corrections mentioned above, there are Higgs boson propagator corrections 
and vertex corrections due to SUSY particles. The first class of radiative corrections, can 
be taken into account by introducing the effective mixing angle a e g that diagonalizes the 
neutral Higgs boson mass matrix [7]. Such type of corrections are known analytically up 
to two-loop order in supersymmetric QCD (SQCD) and the supersymmetric electroweak 
theory (SEW), see [8] and references therein. The second class of radiative corrections are 
especially important for large values of tan/3. Usually, they are derived using the effective 
Lagrangian approach [9]. The one- loop contributions are known since long [6j[T0], while 
the two-loop corrections have been computed very recently jllj . 

It is the aim of this paper to present the complete two-loop SQCD corrections to the 
decay width h — ¥ bb, taking into account the exact dependence on the supersymmetric 
particle masses and working in the full theory. 

The paper is organized as follows: in the next section we introduce our framework and 
the quantities required for the computation of the decay width T(h — > qq) through two 
loops in the MSSM. In Section 3 we describe the actual two-loop calculation, pointing out 
the connection between the radiative corrections to the vertex hqq and those to the quark 
propagator through the low-energy theorem. The formalism discussed in this section 
is valid for a general quark flavour. However, we specify our calculation for the case 
of bottom quark which generates the dominant decay mode of a light Higgs boson. In 
Section 4 we perform a numerical analysis and discuss the phenomenological implications 
of the two-loop SQCD vertex corrections to T(h — > bb). We present the conclusions in 
Section 5. 



^Tor the definition of the parameter xt see Section I27T 



Notation and theoretical framework 



The part of the MSSM Lagrangian describing the fermionic Higgs decay can be written 
in the following form 



£■ — £q C D + ^SQCD + 2_^ £q<t>i + 2_^ £q<t>i 
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where 



£q* = - J2 -zrapQ^i and 4i* = ~J2 J2 -^r^kr& 



(2) 



g =l 



q=l r,k=l,2 



£qcd + £sqcd denotes the supersymmetric extension of the full QCD Lagrangian with six 
quark flavours. The couplings gp and g^ kr are defined in Table [Q m q denotes the mass of 
quark q, v — \fv\ + v\ with Vi,i — 1, 2, the vacuum expectation values of the two Higgs 



doublets of the MSSM. The fields & 



1,2, denote the squark mass eigenstates, while 



q stands for the mixing angle defined through: 
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cot f3 , for up-type quarks 
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where A q is the trilinear coupling and /isusy the Higgs-Higgsino bilinear coupling. The 
fields </>i,i = 1,2, denote the neutral CP-even components of the MSSM Higgs doublets 
and they are related to the Higgs mass eigenstates through the orthogonal transformation 
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As usual, h stands for the lightest Higgs boson. The mixing angle a is determined at the 
leading order through 



tan 2a = tan 2/3 



M\ + M\ 
M\ - Ml ; 
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— < a < 

2 



(5) 



where Mz is the mass of the Z boson and tan/3 = v<lJV\. 

In the following, we assume the mass of the lightest Higgs boson h to be much smaller 
than the mass of the top-quark and of the SUSY particles, as well as all the other Higgs 
bosons. In this case, the physical phenomena at low energies can be described with an 
effective theory containing five quark flavours and the light Higgs. At leading order in 
the heavy masses, the effective Lagrangian Cy can be written as a linear combination of 
three physical operators [12l[T3] constructed from the light degrees of freedom 
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Table 1: Yukawa coupling coefficients for up and down type quark and squark, where 
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where £qcd denotes the Lagrangian of QCD with five active flavours and the coefficient 
functions Cj ,i = 1, 2q, 3q, parametrize the effects of the heavy particles on the low-energy 
phenomena. The superscript labels bare quantities. The three operators are defined as 



o° -- 


/s~i0,',a\2 

- \S*w > > 


W 2q ~ 


= mJW, 


U 3q ~ 


= f>'(ip >'-m° q >')q >', 



(7) 

where G°''£" and D9j' are the gluon field strength tensor and the covariant derivative, 
respectively, and the primes label the quantities in the effective theory. The operator 
Ozq vanishes by the fermionic equation of motion and it will not contribute to physical 
observables. So, the last term in Eq. (jSJ) might be omitted, once the coefficients C°, CL 
are determined. The coefficient functions contain information about the heavy particles 
that were integrated out in the construction of the effective theory. On the contrary, 
as can be understood from Eq. (j7|), the operators encounter only the effects of the light 
degrees of freedom. 

The relations between the parameters and fields in the full and effective theories are 
given by 
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where g s = \Z4ira s is the strong coupling. The coefficients Q"' , Q"' 5 Cg" J , Cm' are the bare 
decoupling coefficients. They may be computed from the transverse part of the gluon 



polarization function and the vector and scalar part of the quark self-energy via [13] 

Cf = i + n°», 
Cf = i + sr (o) , 
C ( 0) = i - ^ ft (o) (9) 



For the derivation of the coefficient Q one has to consider in addition one vertex in- 
volving the strong coupling, for example qqg or ccg, where c denotes the Faddeev- Popov 
ghost. The decoupling coefficients are independent of the momentum transfer, so that 
they can be evaluated at vanishing external momenta. The superscript h indicates that 
in the framework of Dimensional Regularization (DREG) or Dimensional Reduction 
(DRED) only diagrams containing at least one heavy particle inside the loops contribute 
and that only the hard regions in the asymptotic expansion of the diagrams are taken 
into account. They have been computed in QCD including corrections up to the four- loop 
order for the strong coupling [14j and three-loop order for quark masses [13J. In the 
MSSM the two-loop SQCD [T5T - [T7] and SEW [18] expressions are known. Similar to 
the case of SM, the decoupling coefficients derived within the MSSM can be connected 
through the Low Energy Theorem (LET) [TH] with the coefficients C®, CJL. We discuss 

in more detail the relation between the coefficients C^, C° g and £m , Q in Subsection 13.21 

The renormalization procedure of the dimension four operators in the Minimal Sub- 
traction Scheme within DREG (MS) [20] is known since long time p2]. The main aspect 
is that different operators in general mix under renormalization. For the convenience of 
the reader we reproduce the results for the renormalization constants of the operators O® 
and 0\ that are of interest for the fermionic Higgs decays 

Oi = Z u O° + Z 12 0° 2q , C> 2 = Z 22 C>», , where 

z^fi-imy, Za= .^^( 1 .i T §w\-\ Za=l , (10) 



c x = z^cl , C 2q = C° 2q - ^C? . (11) 

The explicit expressions for the /3-function and quark mass anomalous dimension 7 m of 
QCD with ni = 5 active flavours at the one-loop order that are needed in the present 
paper, are given by 
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The bare coefficient functions Cf ,i = l,2q, must be computed diagrammatically For 
the calculation of the 0(a 2 s ) corrections to the process h — > qq, the knowledge of the 
coefficient functions C° and C 2q is required at the one- and two-loop order, respectively. 



The renormalized coefficient functions and operators are finite but not renormalization 
group (RG) invariant. In Ref. [5], a redefinition of the coefficient functions and operators 
was introduced so that they are separately RG invariant. This procedure allows us to 
choose independent renormalization scales for coefficient functions and operators. In 
practice, one makes a separation of scales: one chooses \x ~ Mh for the renormalization 
scale of the operators and /x m M (where M denotes an averaged mass for the heavy 
supersymmetric particles) for the coefficient functions. The new coefficient functions 
read [5] 

C l{ M,M h) = <&^^C l{ M), 
a> s (M h )(3(V(a' s (M)) 



C 2 {M,M, 



h 



- i?; ( ff]>» h£>W&)) - 7gK(M h ))]Ci(M) + C 2q (M) , (13) 
npw{a' a {M)) 

where the superscript (5) marks that n\ = 5 in Eq. (1121) . We employ this approach for the 
evaluation of the decay width F(h — > qq). More details about the practical calculation 
are discussed in the next section. 

2.1 Higgs decay width 

Once the renormalized coefficient functions C\,Ci are known, the decay width for the 
process h — >• qq can be predicted. From Eqs. (jfH) and (GO) one can derive a general formula 
for the inclusive h — > qq decay width [5J 
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where 1^°) represents the complete leading order (LO) result given by 
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As is well known, the large logarithms of the type ln(M^/ra 2 ) can be resummed by taking 

m q in Eq. ( 1T5|) to be the MS mass mj 1 (//) evaluated at the scale \i = Mh- The QCD 
correction A^ CD is known since long time [2~T] . 
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with £(x) being the Riemann's zeta function. The additional QCD correction generated 
through double-triangle topologies S^ CD was first computed in Ref. [5], 

h qcd = 2^0 c ,(_ 1 » + aC( 2)-] 11 «g-61a^) . (17) 



The universal corrections 5 U of 0(a™x t ), where x t = {at/An) 2 = GfM 2 /(8it 2 \/2), with a t 
the top- Yukawa coupling, contain the contributions from the renormalization of the Higgs 
wave function and the vacuum expectation value 
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The coefficient functions C\ ,C 2q (and implicitly C^C-z) are known within SQCD at the 
one-loop order since quite some time [TQ1I23]. For completeness, we display them here 
providing also 0(e) terms that are necessary for the two-loop calculation. 
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where Cp = 4/3, L x = ln(/i 2 /m 2 ) and the functions F\ and F 2 are defined through 
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In the above formulas, a' s (n) denotes the strong coupling constant computed in the MS 
scheme and taking into account n\ = 5 active quark flavours. 

The computation of the coefficient function C 2q through two loops in SQCD is discussed 
in some detail in the next section. 



3 Calculation of the coefficient function Ci q at next- 
to-next-to-leading order (NNLO) 

For the derivation of the coefficient functions one has to compute Green functions in 
the full and effective theory and make use of the decoupling relations Eq. (jSJ) to connect 
them [13]. In general, one Green function contains several coefficient functions. For 
example, the amputated Green function involving the qq pair and the zero-momentum 
insertion of the operator Oh which mediates the couplings to the light Higgs boson h 
contains both coefficient functions C 2q and C$ q . Similarly, one possibility to compute the 
coefficient function C\ involves the Green function formed by the coupling of the operators 
Oh to two gluons. 

In the following, we restrict the discussion to the computation of the coefficient func- 
tion C 2q . Considering the appropriate one-particle-irreducible (1PI) Green function, we 
get 

r° g o>, ~P) =* 2 J dxdye^-^(Tq°(x) ( f(y)O h {0)) 1P1 , (22) 

where p is the outgoing momentum of q. In a next step we express the operator Oh with 
the help of Eqs. (|SJ) and (JTj) and make use of the decoupling relations Eq. (jSJ). One can 
easily see that the above Green function will get contributions only from the operators 
2q and G 3q 

TJ4>, -p) = -Cf J Mye^-y\Tq>'\x)q'»{y){C 2q 2q + C 3q O Sq )) lpl , 

= CfC^^-^V^ + Cf^V- (23) 

In the last step we have used the Feynman rules for the scalar dimension four operators 
that can be found in Ref . [21] and the fact that r=' (p, —p) denotes an amputated 
Green function. Exploiting the fact that the coefficient functions do not depend on the 
momentum transfer, one can set also p = 0. In this case, on the l.h.s. of Eq. (123]) only 
the hard parts of the Green function survive, as the massless tadpoles are set to zero in 
DRED and DREG. As before, the superscript h stands for hard contributions. 

The validity of the approximation m\ = pi ~ was extensively studied in the context 
of the SM and reconfirmed for the case of gluon fusion at two-loop order in SQCD in 
Ref. [25j. We expect that this approximation holds also in the case of fermionic Higgs 
decays, due to the heavy supersymmetric mass spectrum. 

There are two possibilities currently used in the literature for the derivation of 
r=' G (0,0). The first one is the direct computation of the scalar and vector components 
of the vertex function making use of the appropriate projectors. The second one uses the 
LET, which relates the vertex corrections to the hqq coupling to the quark self-energy 
corrections via appropriate derivatives [T§||21)]. 



3.1 Direct calculation of the coefficient function C\ at NNLO 

Decomposing the Green function ^ q : q0h into its scalar and vector components rJ , 
r,' c , v one derives from Eq. (1251) two linearly independent relations for the coefficients 
C% and C 3 ° g 

A \,(o,o) = ^^{cl-cl 



r°^(o,o) = c 2 (0) c 3 V (24) 

In SQCD at the two- loop order there is also an axial contribution to r=' c (0,0), which 
arises from diagrams where a top quark-squark pair is exchanged from a gluino propagator. 
However, it generates only contributions 0(oQ) to T(h — > gg) which are beyond the 
precision we are interested in this paper. 

Finally, using Eqs. © the expression for CL through two loops reads 

r o _ r^ h;a (o,o) r^(o,o) 

29 " 1 - S^(0) + 1 + E&*(0) ' l j 

One can either work in the (0i, 02) basis, i.e. generating two sets of vertex corrections for 
each Higgs field, or one derives the Feynman rules for the couplings of quarks and squarks 
to the light Higgs using the relation Eq. (@J. 

For the computation of the vertex corrections up to two loops we have implemented two 
independent setups. In one of them, the Feynman diagrams are generate with the help of 
the program FeynArts [26], then its output is handled in a self-written Mathematica code 
which includes the two-loop tensor reduction and the mapping of the vertex topologies 
to the two-loop tadpole ones [27J. The last step is possible due to the fact, that we 
neglect the mass of the light Higgs boson and of the external light quarks. In the second 
setup, the Feynman diagrams are generated with the program QGRAF [28], and further 
processed with q2e and exp [29j|30]. The reduction of various vacuum integrals to the 
master integral was performed by a self-written FORM [31J routine. 

3.2 LET derivation of the coefficient function C£, at NNLO 

The connection between the coefficient functions C°, C^ q and the decoupling coefficients 
Csj Cm or equivalently n°' h (0), £°' h (0) and E°' h (0) was extensively studied in the context of 
the SM. The validity of the LET was verified up to three- loop order in QCD [13]. In the 
framework of the MSSM, however, the derivation of LET is much more involved due to the 
presence of the two Higgs doublets and of many massive particles and mixing angles. The 
applicability of LET in SQCD at two-loop order was verified only very recently. Namely, 
the relationship between the coefficient function C\ and the hard part of the transverse 
gluon polarization function n 0,/l (0) has been established in Ref. [25]. Furthermore, the 
leading two-loop contributions to the effective bottom Yukawa couplings have been derived 
from the scalar part of the bottom quark self-energy in Ref. [TT]. It is one of the aims 
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of this paper to verify the relationship between the coefficient function Ciq and the hard 
part of the scalar and vector contributions to the quark self-energy. 

For our calculation it is very convenient to work in (0i,02) basis, which means that 
we have to decompose the Green functions according to Eq. (jlj) 
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Similarly, the coefficient function can be written as follows 



(26) 



U 2q 



sinaC^ +cosaC^ i02 . 



Applying the LETa to the individual components Tggo^-a we get 



(27) 



ggCV ; s 



:s (o,o) 






,(0,0) 



l__d_ 
nib dcf>i 
_d_ 
d<pi 



m 6 (l-^ "(())) 

sr(o) n 



1 



-D n 



<t>i=Vi 



D n 



4>i= v i 



m b 
Eg-*(0) 



m 6 (l - S^(0)) 



(2f 



with i — 1, 2 and a = t>, s. As we are considering only the hard parts of the above Green 
functions no complication related to the occurrence of infrared divergences is encountered. 
In practice, it is convenient to express the operators Dq,^ introduced in Eq. ( )28l) in terms 
of derivatives w.r.t. masses and mixing angles. This can be achieved using the field- 
dependent definition of the parameters, in our case quark and squark masses and squark 
mixing angles [33] . 

The formulas derived up to now are valid for a generic light quark flavour q. However, 
for phenomenological applications the decay channel h — > bb is the most important one. 
The explicit expressions for the operators Db,^ can be easily derived from the Eqs. (11) 
and (12) in Ref. [25\. We quote them here for completeness and to fix our normalization: 
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(29) 



In the above formulas we keep only the terms that do not vanish in the limit mb — > 0. 



2 The gauge-fixing condition for SQCD is independent of the vacuum expectation values i>i,2, so that 
LET holds in its trivial form 132 . 
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As is well known, in Eqs. ([28]) one has first to apply the derivative operators Dq^ and 
afterwards perform the renormalization. For simplicity of the notation we suppress the 
superscript (0), labeling bare quantities. We checked explicitly at the diagram level that 
Eqs. ( 128]) hold through two loops. The computation of the two-loop diagrams contributing 



to Y^ h {Q) goes along the same line as that for T^' q0 ;a (0,0). The exact results together 
with few expansions for special mass hierarchies can be found in Ref. pT6|[TT] . 

Let us mention at this point that for large values of tan (3 the dominant contribution 
to the coefficient function C2 q is contained in the first term in Eq. ( 125]) . The /isusYtan/3- 
enhanced contributions are implicitly resummed in Eq. ( 125]) . through the presence of 
the denominator 1 — X° ,h (0), which contains contribution of the form a™/isusY tan (3. In 
the framework of LET the /isusy tan /9-enhanced contributions to Tg q0h]S are generated 
through the term proportional to the derivative J-j, in Db,cj> 2 - Taking into account the 
parametric dependence of E°' on masses and mixing angles, one can easily derive these 
contributions from the terms proportional to sin 26^ in S^Jj Such contributions have also 
been derived in Ref. [IT] using the effective Lagrangian approach. Indeed, after discarding 
the additional pieces comprised in our computation of the coefficient C^, namely the 
vector part and the rest of the derivative operators in Eq. ([29]) . we get good numerical 
agreement with the results of Ref. [IT] . 



3.3 Regular izat ion and renormalization scheme 

It is well known that the appropriate regularization scheme for the computation of ra- 
diative corrections in supersymmetric theories is DRED. However, the most convenient 
regularization scheme for the handling of the dimension four operators at higher orders 
is DREG as discussed in Section [2j For the renormalization we employed two differ- 
ent approaches. In one of them we computed the radiative corrections to the coefficient 
functions directly in DREG. This implies that some of the supersymmetric relations be- 
tween couplings of quarks and squarks do not hold anymore. More precisely, one has 
to distinguish between the gluino-quark-squark coupling g s and the gauge coupling g s , 
and between the Yukawa couplings of Higgs bosons to quarks gp and squarks gp kr - The 
relationships between the different couplings are necessary only at the one-loop order and 
they are well known since long time [51] . 

In the second approach, we performed the two-loop computation of the coefficient func- 
tions in the DRED scheme and afterwards converted the results into the DREG scheme 
using the two-loop translation relations for the quark masses and strong couplings de- 
fined in the full [35] and effective theory [36], and Eqs. ( )28l) . As a consistency check, we 
explicitly verified that the results obtained with the two methods agree. 

For the renormalization of the divergent parameters we used the on-shell scheme for the 
gluino and bottom squark masses and mixing angle and the minimal subtraction scheme 
MS or DR for the strong coupling and the bottom quark mass. The renormalization of the 
trilinear coupling A^ was performed implicitly through the use of relation Eq. (|3]). The 



5 A similar observation was made in Ref. [II] . too. 
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explicit formulas for the one-loop counterterms are well-known in the literature (see for 
example Ref. [37]). 

The complete two-loop results for the SQCD corrections to the coefficient function 
C*2g discussed in this section are too lengthy to be given here. They are available in 
MATHEMATICA format upon request from the authors. For further applications, we also 
provide results for the case where all parameters are renormalized minimally. 

In principle, the results obtained in this section can be easily generalized to the heavy 
Higgs decays taking into account the necessary changes in the Yukawa couplings as given 
in Table [U However, the application of effective theory formalism introduced in Section [3] 
is not justified in this case, i.e. the condition M H <C M t , M SU sy does not hold anymore. 

4 Numerical results 

In this section we study the phenomenological implications of the two-loop corrections to 
the coefficient function Ci q on the Higgs decay width T(h—> bb). The SM input parameters 
are the strong coupling constant at the Z-boson mass scale a s (Mz) = 0.1184 [38J, the 
top quark pole mass M t = 173.1 GeV [39] and the running bottom quark mass in the MS 
scheme m^rrif,) = 4.163 GeV [40]. For the supersymmetric mass spectrum we adopted 
the corresponding values of the "small a e ff" and "gluophobic" scenarios as defined in 
Ref. 00 

For the running of the strong coupling constant within QCD we use the Mathematica 
package RunDec [41]. For the evaluation of the strong coupling constant within the six- 
flavour SQCD and the DR scheme we follow Ref. [PT]rl 

An important ingredient for the computation of the decay width T(h — > bb) is the effective 
mixing angle of the neutral Higgs sector a eff [7] that takes into account the radiative 
corrections to the Higgs propagator. In practical applications one replaces the tree-level 
mixing angle a defined in Eq. (J5j) with a e g. This can be computed in perturbation theory 
from the knowledge of the radiative corrections to the self-energy matrix of the neutral 
Higgs doublet E^, S 02 , S 01</>2 , 

= -(Ml + MDsin/^-^ (30) 

Mf cos 2 (3 + M\ sin 2 (3 - M 2 - S 01 

where M^ stands for the on-shell mass of the light Higgs boson. For the numerical analyses 
we implemented the exact two-loop results from Ref. [42] rl 

In Fig. [1] we show separately the renormalization scale dependence of the coefficient 
functions Ciq^ and C2q,<p 2 , which describe the effective Yukawa couplings of the neutral 
Higgs fields 4>i an d 02- We choose the "small a e ff" scenario with tan/3 = 50, Ma = 



4 We used the tree-level formulas to derive the mass eigenvalues for squark fields. 
5 For a more detailed discussion see Ref. [17] and the references therein. 

6 Very recently the three-loop SQCD corrections to Mh have been computed [43]. However, they are 
valid only for specific mass hierarchies. 
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Figure 1: The renormalization scale dependence of the coefficient functions 
C2q,(j>-L cos (3 (a) and Ci q ,fa sin (3 (b) is depicted at one- and two-loop order in the 
"small a e ff" scenario. 
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300 GeV and evaluate a e g with the tree-level formula in order to avoid additional scale 
dependence. The dashed and solid lines correspond to the one- and two-loop results, 
respectively. As can be read from the figure, at the one-loop order the scale dependence 
amounts to about 50% when the renormalization scale is varied around the average value 
/io = (rag + m~ b2 + ra§)/3 ~ 658 GeV by a factor 10. At the two-loop order the variation 
with the renormalization scale is significantly improved. The remaining scale dependence 
is below 6%. An interesting aspect is the opposite evolution of the two coefficient functions 
with the renormalization scale. This feature is reflected in a milder scale dependence 
of the full coefficient function C 2q of about 4% and 0.5% at one- and two-loop order, 
respectively. However, for low A-boson masses Ma < 120 GeV this numbers change to 
27% and 5%, respectively. As usual, we can interpret the last number as an estimation of 
the theoretical uncertainties due to unknown higher order corrections. So, the two-loop 
SQCD corrections are essential for the accurate prediction of the decay width F(h—> bb). 

In Fig. [2] we display the renormalization scale dependence for the "gluophobic" sce- 
nario, where we fixed the value of tan/3 to 20, Ma = 300 GeV and maintain the same 
convention for the lines. A similar behaviour as for the "small a e s" scenario is observed. 
However, in this case the coefficient function C 2q ^ 2 has a much stronger scale dependence 
at the one- loop order than the coefficient C 2q ^ x ■ The scale variation of the full coefficient 
function C 2q sums up to 1.5% and 0.2% for Ma = 300 GeV at the one- and two-loop order, 
respectively. However for Ma < 120 GeV the scale variation amounts to 8% and 1.5% at 
one- and two-loop order, respectively, which shows the importance of the two-loop SQCD 
corrections for this region of the parameter space. 

We can conclude that for phenomenological analyses the choice of the renormalization 
scale around /i ensures small radiative corrections and a good convergence of the pertur- 
bative scale. In the following, we set \i ~ /i for the computation of the SQCD corrections 
to the coefficient functions C 2q and the decay width F(h -> bb). 

In Fig. [3] we depict the dependence of the full coefficient function C 2q on the A boson 
mass Ma (a) and the light Higgs boson mass Mh (b) for the "small a e s" scenario. We set 
tan = 50 and vary the mass of the A boson between 100 GeV < Ma < 200 GeV. For the 
evaluation of the Higgs boson mass and the effective mixing angle a e ff we employed the 
two-loop SQCD results [42J. As can be seen from Fig. E^a) for low Ma values there are 
large one-loop corrections of about 60% of the tree-level values. They originate from the 
large corrections to the scalar part of the quark self-energy, that are actually resummed 
through the use of the formula given in the Eq. (j25l) . This feature is also reflected by 
the relatively small two-loop corrections of about 2% of the tree-level values. For Ma 
values larger than 130 GeV one observes a steep increase of the coefficient function C 2q 
(decrease of the absolute value) due to cancellation of Eq. (130]) . that implies a c s — ¥ 0. In 
this case, C 2q reaches its minimal absolute value. From panel (b) one notices a similar 
steep increase of the coefficient function C2 q when Mh reaches its maximal value of about 
125 GeV. A similar behaviour is also observed for the "gluophobic" scenario. So, for both 
scenarios we expect a strong suppression of the decay width T(h — > bb) for large A-boson 
masses or equivalently for M^ close to its maximal value for which a e fj — > 0. 

In Fig. H] the dependence of the coefficient function C 2q on tan (3 is shown for the 
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Figure 2: The renormalization scale dependence of the coefficient functions 
C2q,(j>i cos j5 (a) and Ci q ,fa sin (3 (b) is depicted at one- and two-loop order in the 
"gluophobic" scenario. 
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"small a e ff" (a) and "gluophobic" (b) scenarios, where the A-boson mass was fixed to 
Ma = 130 GeV. As expected, we observe a significant increase of the magnitude of the 
radiative corrections with the increase of the tan/3 value. At the one-loop order the 
radiative corrections amount to about 70% (a) and to 33% (b) from the tree level values. 
At the two-loop order they sum up to 3% and 5%, respectively. 

In Fig. [5] we display the decay width for h — > bb as a function of the Higgs boson mass 
Mh, considering the "small a e ff"(a) and "gluophobic" (b) scenarios. We chose tan/3 = 50 
and tan fl = 20 for the case (a) and (b), respectively. The two- loop genuine QCD and 
EW corrections to the process h — > bb, as well as the two-loop SQCD corrections to the 
Higgs boson propagator are depicted by the dotted lines. More precisely, they are derived 
from Eq. (Ib£|) . where the coefficient functions CiandC 2 are set to their tree-level values. 
The additional SQCD vertex corrections parametrized through the coefficient functions 
C\ and C 2 are represented at the one- and two-loop order by the dashed and solid lines, 
respectively. We also take into account the one-loop SEW corrections to the coefficient 
function C 2 and fix their renormalization scale at //sew — ( m t! + m i 2 + / i susY)/15, for 
which the two-loop SEW corrections become negligible [TT] . 

For a relatively light Higgs boson mass Mh-, the large one- loop radiative corrections of 
about 70% (a) and 50% (b) are still amplified by mild two-loop corrections that can 
reach as much as about 8% from the decay width including QCD corrections even for 
the selected choice of the renormalization scale of SQCD corrections. The large SQCD 
radiative corrections to T(/i -)• bb) have only a relatively small impact on the branching 
ratio BR(h — > bb) but they can have a large impact on BR(h —¥ t + t—). For sufficiently 
large tan/3 and /isusY, the measurement of BR(h — > t + t—) can provide information 
about the distinction between the SM and MSSM predictions. 

For a large Higgs boson mass for which a e Q — > 0, the partial decay widths for h — >■ bb 
and h — > t + t~ are significantly suppressed. In this case the radiative corrections (in 
particular the corrections to the Higgs boson propagator in Eq. (130]) ) are essential for an 
accurate prediction of T(h — > bb) and T(h — > t + t~). Furthermore, the BR(h — > 77) will 
be strongly enhanced, improving the LHC prospects of finding a light Higgs. 

5 Conclusions 

The knowledge of the Higgs boson couplings is essential for its searches at the present 
hadron colliders. In this paper we calculate the two-loop corrections of 0(a 2 s ) to the 
Yukawa couplings within the MSSM. We employed the effective Lagrangian approach un- 
der the assumption of large top quark and supersymmetric particle masses. We calculate 
analytically the two-loop corrections to the coefficient function C 2(? , taking into account 
the complete mass dependence. For large values of tan /3 the radiative corrections need to 
be resummed, which in our approach is performed through the use of the formula given in 
Eq. (|2"5]) . Furthermore, we verified at the diagram level the applicability of the low-energy 
theorem for Higgs interactions in the framework of the MSSM as stated in Ref. [2"5] . 
From the phenomenological point of view, the two-loop corrections presented here 
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reduce significantly the theoretical uncertainties, estimated through the variation with 
the renormalization scale, at the percent level. The two-loop SQCD corrections become 
sizable for tan/3 > 20 and M A < 130 GeV. 
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